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Abstract 

Motivated by the recent observation of CP violation in the charm sector by LHCb, we study 



direct CP asymmetries in the standard model (SM) for the singly Cabibbo-suppressed two-body 
hadronic decays of charmed mesons using the topological-diagram approach. In this approach, 
the magnitude and the phase of topological weak annihilation amplitudes which arise mainly from 
final-state rescattering can be extracted from the data. Consequently, direct CP asymmetry ajj^. 
at tree level can be reliably estimated. In general, it lies in the range 10~ 4 < a!^^ < 10 -3 . 
Short-distance QCD penguins and penguin annihilation are calculated using QCD factorization. 
Their effects are generally small, especially for D — > VP modes. Since weak penguin annihilation 
receives long-distance contributions from the color-allowed tree amplitude followed by final-state 
rescattering, it is expected to give the dominant contribution to the direct CP violation in the 
decays D° — > K + K~ and D° — > ir + ir~ in which is absent. The maximal Aajfep, the direct 

CP asymmetry difference between the above-mentioned two modes, allowed in the SM is around 
—0.25%, more than 2a away from the current world average of —(0.645 ± 0.180)%. 



1 



I. INTRODUCTION 



Recently the LHCb Collaboration has reported a result of a nonzero value for the difference 
between the time-integrated CP asymmetries of the decays D° — > K + K~ and D° — > 7t + tt~ [l| 

AA C P = Acp(K + K-)-A C p{7r + 7T~) = -(0.82 ±0.21 ±0.11)% (LHCb) (1) 

based on 0.62 fb" 1 of 2011 data. The significance of the measured deviation from zero is 3.5a. 
However, based on a data sample corresponding to the integrated luminosity of 5.9 fb _1 , the CDF 
Collaboration 0] obtained Acp(tt + tt-) = (0.22 ± 0.24 ± 0.11)% and A C p(K+K-) = -(0.24 ± 
0.22 ± 0.09)%, and hence 

AA CP = -(0.46 ±0.31 ±0.11)% (CDF). (2) 

The time-integrated asymmetry can be written to first order as 

AcpU) = a^ P (f) + { ^Sp , (3) 

where a^p is the direct CP asymmetry, a^p is the indirect CP asymmetry, (t) is the average 
decay time in the sample and r is the D° lifetime. A global fit to the data of AAcp is consistent 
with no CP violation only at 0.128% CL Q]. The central values and ±1<7 errors for the individual 
parameters are: a£ P = -(0.019 ± 0.232)% and Aa d j P = -(0.645 ± 0.180)%. 

Whether the first evidence of CP violation in the charm sector observed by LHCb is consistent 
with the standard model (SM) or implies new physics will require further analysis of more data and 
improved theoretical understanding. For some early and recent theoretical studies, see Refs. [j- 14]. 

It is known that in order to induce direct CP violation, one needs at least two distinct decay 
amplitudes with non-trivial strong and weak phase differences. In B physics, there exist several 
QCD-inspired approaches describing the nonleptonic B decays, such as QCD factorization (QCDF) 



151 ]. pQCD [16[ and soft collinear effective theory [17[. Even so, one needs to consider 1/mj power 



corrections in order to explain the observed direct CP asymmetries in B decays. For example, in 
QCDF it is necessary to take into account two different types of power correction effects in order 
to resolve the CP puzzles and rate deficit problems with penguin-dominated two-body decays of B 
mesons and color-suppressed tree-dominated 7r°7r° and p°TT° modes: penguin annihilation and soft 
corrections to the color-suppressed tree amplitude [18]. 

The situation is far worse in the charm sector as a theoretical description of the underlying 
mechanism for exclusive hadronic D decays based on QCD is still not yet available. This is 
because the mass of the charmed quark, being of order 1.5 GeV, is not heavy enough to allow for a 
sensible heavy quark expansion. Indeed, it does not make too much sense to generalize the QCDF 
and pQCD approaches to charm decays as the l/m c power corrections are so large that the heavy 
quark expansion is beyond control. In short, there is no reliable model which allows us to estimate 
the phases and magnitudes of the decay amplitudes beyond the color-allowed tree amplitude. 

Luckily, we do have a powerful tool which provides a model-independent analysis of the charmed 
meson decays based on symmetry, namely, the diagrammatic approach. In this approach, the 
flavor-flow diagrams are classified according to the topologies of weak interactions with all strong 
interaction effects included. Based on flavor SU(3) symmetry, this model-independent analysis 
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enables us to extract the topological amplitudes and probe the relative importance of different 
underlying decay mechanisms. It is complementary to the factorization approaches. Analysis 
based on the flavor-diagram approach indicates a sizable weak annihilation (^-exchange or W- 
annihilation) topological amplitude with a large strong phase relative to the tree amplitude. Since 
weak annihilation and final-state interactions (FSI's) are both of order l/m c in the heavy quark 
limit, this means FSI's could play an essential role in charm decays. Indeed, weak annihilation 
contributions arise mainly from final-state rescattering, and this explains why an approach based 
on heavy quark expansion in l/m c is not suitable for charm decays. 

The great merit and the strong point of the topological-diagram approach is that the magnitude 
and the relative strong phase of each individual topological tree amplitude in charm decays can be 
extracted from the data. 1 This allows us to calculate CP asymmetry at tree level in a reliable 
way, granting us an idea about the size of CP violation in charmed meson decays. 

Since there is no tree-level direct CP violation in D — > K + K~ and D — > tt + tt~ decays, 
the observation of AAqp by LHCb indicates that, contrary to the conventional wisdom, penguin 
diagrams in singly Cabibbo-suppressed (SCS) decay channels do play a crucial role for CP violation 
even though they may not affect the branching fractions. Indeed, this observation is quite natural 
in the topological approach since, just as the enhancement in weak annihilation diagrams through 
FSI's, weak penguin annihilation (more specifically, the QCD-penguin exchange diagram to be 
introduced in the next section) also receives contributions from the color-allowed tree amplitude 
followed by final-state rescattering. 

The time- integrated CP asymmetry receives indirect CP -violating contributions denoted by a™p 
through the D°-D mixing. Such a mixing is governed by the mixing parameters x = (m\ — m-zj/T 
and y = (Ti — r2)/(2r) for the mass eigenstates D\ and D2, where mi^ and T\p are respective 
masses and decay widths and T is the average decay width. It is known that the short-distance 



contribution to the mixing parameters is very small 19], of order 10 6 . On the theoretical side, 



there are two approaches: the inclusive one relying on the l/m c expansion (see e.g ., Ref. 20J] for a 



recent study), and the exclusive one with all final states summed over. In Ref. 2l|], only the SU(3) 
breaking effect in phase space was considered for the estimate of y. Consequently, the previous 
estimate of mixing parameters is subject to larger uncertainties. We believe that a better approach 
is to concentrate on two-body decays and rely more on data than on theory. This is because the 
measured two-body decays account for about 75% of hadronic rates of the D mesons. For PP and 
VP modes, data with good precision for Cabibbo-favored (CF) and SCS decays are now available. 
For as-yet unmeasured doubly Cabibbo-suppressed (DCS) modes, their rates can be determined 
from the diagrammatic approach. We obtained x = (0.10 ± 0.02)% and y = (0.36 ± 0.26)% from 
the PP and VP final states which account for nearly half of the hadronic width of D° 22|] . It is 



conceivable that when all hadronic states are summed over, one could have x ~ (0.2 — 0.4)% and 
y ~ (0.5 — 0.7)%. At any rate, indirect CP violation is suppressed by the smallness of the mixing 
parameters x and y and, moreover, it tends to cancel in the difference between K + K~ and 7r + 7r _ 



1 This is not the case in B decays where one has to make additional assumptions in order to extract 
individual topological amplitude cleanly from the data. 
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final states. 

There are also a lot of experimental efforts measuring direct CP asymmetries in D decays with 
3-body and 4-body final states. This will involve a more complicated and time-consuming Dalitz- 
plot analysis. Nevertheless, the analysis will be very rewarding as it will provide us with much 
more information about the underlying mechanism for CP violation. We plan to explore this topic 
elsewhere. 

The purpose of this work is to provide a realistic SM estimate of direct CP violation in two- 
body hadronic D decays based on the topological-diagram approach. In Sec. II we recapitulate the 
essence of the diagrammatic approach and its application to SCS D — > PP and D — > VP decays. 
The importance of final-state interactions is emphasized. The topological amplitudes are related 
to the quantities in the QCDF approach in Sec. III. Sec. IV is devoted to the calculation of direct 
CP asymmetries by taking into account SU(3) breaking and penguin effects. Sec. V contains our 
conclusions. 



II. DIAGRAMMATIC APPROACH 
A. Analysis for hadronic charm decays 

It has been established sometime ago that a least model-dependent analysis of heavy meson 



decays can be carried out in the so-called quark-diagram approach 2314251] . In this diagrammatic 



scenario, the topological diagrams can be classified into three distinct groups as follows (see Fig. [T]): 

1. tree and penguin amplitudes: 

• T, color-allowed external VF-emission tree amplitude; 

• C, color-suppressed internal ^-emission tree amplitude; 

• P, QCD-penguin amplitude; 

• S, singlet QCD-penguin amplitude involving SU(3)F-singlet mesons (e.g., 77W, u>, 4>); 

• Pew 5 color- favored EW-penguin amplitude; and 

• P^w, color-suppressed EW-penguin amplitude; 

2. weak annihilation amplitudes: 

• E, W-exchange amplitude; 

• A, VF-annihilation amplitude; (E and A are often jointly called "weak annihilation" 
amplitudes.) 

• PE, QCD-penguin exchange amplitude; 

• PA, QCD-penguin annihilation amplitude; 

• PEewj EW-penguin exchange amplitude; and 

• PAeW) EW-penguin annihilation amplitude; (PE and PA are also jointly called "weak 
penguin annihilation" amplitudes) and 
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FIG. 1: Topology of possible flavor diagrams: (a) color-allowed tree T, (b) color-suppressed 
tree C, (c) QCD-penguin P, (d) singlet QCD-penguin S diagrams with 2 (3) gluon lines for M2 
being a pseudoscalar meson P (a vector meson V), where M2 is generally the emitted meson or 
it contains a quark from the weak vertex for the annihilation diagram, (e) VF-exchange E, (f) W- 
annihilation A, (g) QCD-penguin exchange PE, and (h) QCD-penguin annihilation PA diagrams. 
The color-suppressed E W-penguin -P^w an< ^ color-favored EW-penguin Pew diagrams are obtained 
by replacing the gluon line from (c) and all the gluon lines from (d), respectively, by a single Z- 
boson or photon line. The EW-penguin exchange PEew an d EW-penguin annihilation PAew 
diagrams are obtained from (g) and (h), respectively, by replacing the left gluon line by a single 
Z-boson or photon line. The gluon line of (e) and (f) and the right gluon line of (g) and (h) can 
be attached to the quark lines in all possible ways. 



3. flavor-singlet weak annihilation amplitudes: all involving SU(3)F-singlet mesons, 

• SE, singlet ^-exchange amplitude; 

• SA, singlet ^-annihilation amplitude; 

• SPE, singlet QCD-penguin exchange amplitude; 

• SPA, singlet QCD-penguin annihilation amplitude; 

• SPEew, singlet EW-penguin exchange amplitude; and 

• SPAeWj singlet EW-penguin annihilation amplitude. 



261 ] for details. 



The reader is referred to Ref. 

It should be stressed that these diagrams are classified purely according to the topologies of 
weak interactions and flavor flows with all strong interaction effects encoded, and hence they are 
not Feynman graphs. All quark graphs used in this approach are topological and meant to include 
strong interactions to all orders, i.e., gluon lines and quark loops are included implicitly in all 
possible ways. Therefore, analyses of topological graphs can provide information on FSI's. 

The decomposition of the decay amplitudes of SCS D — > PP and D — > VP modes in terms of 
topological diagrams is displayed in Tables [Tl— fTTl Since we will concentrate exclusively on SCS D 
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TABLE I: Branching fractions and invariant amplitudes for singly Cabibbo-suppressed decays of 
charmed mesons to two pseudoscalar mesons. It is understood that the amplitudes with the CKM 
factor X p = V* p V up are summed over p = d,s. Data are taken from 



23]- Predictions based on our 



best-fitted results in @ with exact flavor SU(3) symmetry are given in the last column. 





Mode 


Representation 


«-'exp 

(xW~ 3 ) 


^theory 

fxlO -3 ") 

I 1 


D° 




A \(T 4 E)S ,1 4 P p 4 PE 4 PA] 


1 400 ± 026 


2 24 ± 10 




vr 7r° 


J_ \ \<—C 4 E)5 nr i 4 P p 4- PE 4 PA)} 


80 4- 05 

KJ • KJKJ — 1 — \J • \J t_» 


1 35 ± 05 




'/ 


XJ—E5^ cos 6 — \=C5™ sin 6 4 (P p 4- PE) cos oil 


68 ± 07 


75 ± 02 




7rV 


A P [-P<5 pd sin <p 4 ^Ctfps cos (/) + (P p + PE) sin 0] 


0.89 4 0.14 


0.74 ± 0.02 




'/'/ 


J_ \ S\((l 4- FAft j -A- P p -A- PE 4- PA\ ens 2 rA 
~r{~ ~^2^ sinzo) -f- LUi sin cp)O ps \ 


i «7 _i_ n on 


1 44 4 D 08 

j-.'-t'-t _1_ u.uo 




ijrj' 


4 u cos z<p — hi sm ze? J o ps } 


l.UO =t u.zo 


i in _l n n7 
i.iy =t u.u/ 




A A 


A p ](T 4 Ajdps 4 P p 4 -Pa 4 PA) J 


o.96 ± 0.08 


1.92 ± 0.08 




A A 


Ap(A p 4 2AA) 


U.6az ± U.116 


n 
U 


D + 


7T+7T 


^A d (T4C) 


1.19 4 0.06 


0.88 ±0.10 






a p [-^{1 4 O 4 ZA)dpd cos 03 — UOps sine? 
+V2{P p 4 PP) cos 0] 


Q P^Q _i_ n oi 
o.oo ± U.zl 


1.4o =b U.zo 




ir + r}' 


AJ-WT 4 C 4 2 A) 5 pd sin 4- C5 ps cos ^ 


4.67 4 0.29 


3.70 ±0.37 




K+K° 


+V2(PP 4 PP) sin 0] 








A p LA5 pd 4T5 ps 4P p ±PP] 


5.66 4 0.32 


5.46 ±0.53 


Df 


ir+K° 


A p [T5 pd ±A5 ps 4P p ±PP] 


2.42 ±0.16 


2.73 ±0.26 




7T°K+ 


^A p [-C7,5 pd 4 A5 ps + P p + PE] 


0.62 ±0.21 


0.86 ± 0.09 




K + 7] 


±\plC5 pd 4 A5 ps + P p + PE] cos 

-A P [(T + C + A)5 ps + P p + PE] sin 


1.75 ±0.35 


0.78 ±0.09 




K + T]' 


^A p [C5 pd 4 ^<5 ps +P + P.E] sin 

4A P [(T + C + A)5 ps + P p + PE] cos a6 


1.8 ±0.6 


1.07 ±0.17 



"The subscript p in E' refers to the quark-antiquark pair popping out of the vacuum in the final state. 



decays, the primes on the amplitudes given in Ref. [29(] are dropped. For simplicity, flavor-singlet 
weak annihilation and weak penguin annihilation amplitudes are neglected in these two tables. 

The topological amp litudes T,C, E, A are extracted from the CF D —> PP decays to be (in 
units of 10" 6 GeV) Q (see also H) 

-i(152±l)° 



T 
E 



3.14 ±0.06, 



(1.53i^ 8 V( 122±2 )°, A=(0.39lg:Jgje 



C 
A 



2.61 ±0.08) e" 



i(31^0)c 



(4) 
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TABLE II: Same as Table [J except for singly Cabibbo-suppressed decays of charmed mesons to one 
vector and one pseudoscalar mesons. Due to the lack of information on Ap and Ay, no prediction 
is attempted for D + and decays except D + — > tt + (/). 



Mode 



Representation 



A P [(TV + E P )8 pd ± P£ ± PA P + 
\ P [(T P + E v )5 pd + Pp + PA V + 
2 Cp — Cy + E P + E v )5 pd 

+P P + P v + PAp+~ 
\ p [{T v + Ep)5 pa + pP+^ 
Ap[(T P + E v )S p3 + P P + PE V + PA V ] 
\ p (E v 5 pd + Ep5 ps + PA P + PA V ) 
"1 P ' 



PE P ] 
PE V ] 

PE P - 
PA P ] 



Boxp 


Bti,™,.,, fA,A11 


^theory (^j^lj 


(xicr 3 ) 


(xio- 3 ) 


(xl0~ 3 ) 


4.96 ±0.24 


3.92 ±0.46 


5.18 ± 0.58 


9.8 ±0.4 


8.34 ± 1.69 


8.27 ± 1.79 


■3 79 _|_ n 99 


z.yo in u.yo 


o.o*± nz u.oo 


1.56 ±0.12 


1.99 ±0.24 


1.99 ± 0.22 


4.38 ± 0.21 


4.25 ± 0.86 


3.18 ± 0.69 


< 1.5 


0.29 ±0.22 


0.05 ± 0.06 


< 0.84 


0.29 ±0.22 


0.05 ± 0.06 


< 0.26 


0.10 ±0.18 


1.01 ± 0.18 


1.31 ± 0.17 


1.22 ±0.08 


1.11 ± 0.05 


2.21 ± 0.23 a 


3.08 ± 1.42 


3.94 ± 0.61 




0.07 ±0.02 


0.15 ± 0.01 


0.14 ±0.05 


0.31 ± 0.10 


0.41 ± 0.08 




1.11 ± 0.86 


1.17 ± 0.34 




0.14 ± 0.02 


0.26 ±0.02 







7T 

7r°p° 



K+K 
K~K 
K°K* 
K°K*° 



IjLO 



T] U) 



-PE V ] 



"p\J^V u pd i ^pvps -r J ^ 

\ P {E P 5 pd + E v S ps + PA P + PA V ) 
\\ p l(-C v + C P ~Ep~ E v )8 pd +P p + Pi- ± PE P + PE 



V<t> 
VP 



v'p° 



PA V ) 

Qpd + tiyd ps + fAp + PAy) 
|A p [(— Cy + Cp — Ep — Ey) 

73 AsCp 

^\p[(Cy + Cp + Ey + Ep)S pd cos</> - Cy8 pa sin 4> 

+ (Pp +Py+ PEp + PEy + PA P + PAy) COS </>] 
^\ p [(Cy + C P + Ey + Ep)8 pd s'mrj) + C v 8 ps cose/) 

+ (P p + P V + PEp + PEy + PA P + PAy) sin 0] 
\ p {[-±=Cp cos <j>- (E v + Ep)sin4>]6 pa + (PA P + PAy) sin « 
\\p[{Cy — Cp — Ey — Ep)S pd cost/) — \/2CyS pa s'mcf) 

+ (Pp +Py + PEp + PEy) cos <j>] 
\^p[{Cy - Cp - Ey - E P )S pd sin0 + \/2CyS pa cos</> 
+ (Pl + P v + PEp + PE,A sin<Al 



D+ TT + p" 

n°p+ 



-^\ P [{Ty +C P -Ap + Ay)8 pd + P v -pP+ PEp - PEy] 
^=A p [(Tp +Cy+A P - Ay)8 pd + P p p - Py + PEy - PEp] 

^=A P [(7V +Cp + A p + A v )5 pd + p p+ P v + PEp + PEy] 
X s Cp 

-±=Ap[(T P + C v + Ay + A P )5 pd cos0 - \/2C v 5 ps sin</> 
-(P P + P v + PEp + PEy) cos < " 



0.81 ± 0.15 



7T 1 LJ 



< 0.34 
5.42 



0.16 
0.18 



.21 ±0.43 



V'p + ^A p [(Tp + Cy + Ay + A P )8 pd sin(j) + y/2C v 8 ps cos<f> 



K+K'° 



Df -k+K* 
n K*+ 

K+p° 
K°p+ 
r]K*+ 

r,'K*+ 

K + LU 

K+4> 



-(P P +P V + PEp + PEy) sin < 
MAvS pd + T v S ps + P V + PEp) 
MApS pd + T P 5 pa +P P + PEy) 



3.681^ 
32 ± 11 



5.68 ±0.28 



M T vS pd + AyS pa + P V + PEp) 
-±\p [CyS pd - A v 5 pa -Py- PEp] 
^A p [c P 5 pd - A P 5 pa -P P P - PEy] 
Xp(T P 6 pd + ApS ps + P p + PEy) 

■^Xp{(CyS pd + AyS ps + Py + PEp) cos</> 

-[(Tp + Cy + Ap)S ps +P p + PEy] sin</)} 
^=Ap{ {CyS pd + AySps + Pi- + PEp) sm<j> 

-l(T P +Cy+ A P )S pa +P p + PEy] cos 
-±=Ap [Cp5 pd + Ap8p 3 +P p + PEy] 

Ap[(T v + C P + Ay)8 ps +P V + PEp] 



2.25 ±0.39 



2.7 ± 0.5 



< 2.4 

< 0.6 



"Data from 



2N 



for 



40.4 C 



3l| , where 4> is the r\ — r{ mixing angle defined in the flavor basis 



(5) 



/ T] \ _ / cos <p - sin (p \ I rj q \ 
\r]' J \ sin 4> cos 4> J \ i] s J ' 

with i] q = -^{uu + dd) and 7] s = ss. 

For D — >• VP decays, there exist two different types of topological diagrams since the spectator 
quark of the charmed meson may end up in the pseudoscalar or vector meson. For reduced am- 
plitudes T and C in D — > VP decays, the subscript P (V) implies that the pseudoscalar (vector) 
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meson contains the spectator quark of the charmed meson. For E and A amplitudes with the 
final state q\q~2, the subscript P (V) denotes that the pseudoscalar (vector) meson contains the 
antiquark qi- 

There are two different ways of extracting topological amplitudes: either 

u pot. 

by summing over the the polarization states of the vector meson, or through the relation 

by taking the polarization vector out of the amplitude, where A = (my /mp,)A(e ■ pd)- There 
exist two solutions, denoted by (A) and (S), for the Ty, Cp, and Ep amplitudes, depending on 
whether Eq. ([6]) or Eq. ([7]) is used to extract the invariant amplitudes. Assuming that Tp and Ty 



are relatively real, we obtain two best solutions (Al) and (SI) for Tp, Cy, and Ey 29J (see also 



Ref. 



m 



(A) Ty = 4.lQt° f 7 , Cp = (5.Ut il)e i{162±3) °, E P = (3.09 ± 0.11) e -^ 93±5 )° , 

(S) T V = 2.15i°;°9, C P = (2.68±g;i|) e*( 164±3 )° , E P = (1.69 ± 0.06) e -*( 103±4 )° , 
(Al) T P = 8.11^;1, Cy = (4A5±°f 7 )e i ^ + -*o)°, Ey = (l.Slt^V^ 124 ^) , 
(SI) T P = 3.14±0;|i, Cy = (j.33+0.36) j{vn™)<> ^ Ey = (i.3i+o.40 )e -(io6lf 6 )° (g) 

Solutions (A) and (Al) in units of 10 -6 are obtained using Eq. ([7]), while solutions (S) and (SI) 
in units of 10 _6 (e ■ pd) are extracted using Eq. (|6|). Note that all solutions are extracted from the 
CF D — > VP decays and that the Ap and Ay amplitudes cannot be completely determined based 



on currently available data 29]. 



Under the flavor SU(3) symmetry, one can use the topological amplitudes extracted from the CF 
modes to predict the rates for the SCS and DCS decays. The branching fractions of SCS D — > PP 
decays predicted in this way are shown in the last column of Table HI In general, the agreement 
with experiment is good except for discrepancies in some modes. For example, the predicted rates 
for 7T + 7r _ and 7r°7r° are too large, while those for K + K~ , 7T + r]^ and K + rj^ are too small compared 
to experiments. The decay D° — > K°K° is prohibited by the SU(3) symmetry, but the measured 
rate is comparable to that of D° — > 7r°7r°. 

Some Cabibbo-suppressed modes exhibit sizable violation of flavor SU(3) symmetry. We find 
that part of the SU(3) breaking effects can be accounted for by SU(3) symmetry violation mani- 
fested in the color-allowed and color-suppressed tree amplitudes. However, in other cases such as 
the ratio R = T(D° — > K + K~)/T(D° — > ir + ir~), SU(3) breaking in spectator amplitudes is not 
sufficient to explain the observed value of R. This calls for the consideration of SU(3) violation in 
the PU-exchange amplitudes. 

B. Final-state rescattering 

From Eq. (JJ|) we see that the color-suppressed amplitude C is not only comparable to the tree 
amplitude T in magnitude but also has a large strong phase relative to T. (It is 180° in naive 



S 



factorization.) The VF-exchange E is sizable with a large phase of order 120°. Since ^-exchange 
is of order l/m c in the heavy quark limit, this means that l/m c corrections are very important in 
charm decays. Finally, we see that W- annihilation is substantially smaller than W-exchange and 
almost perpendicular to E. 

In naive factorization, the factorizable weak annihilation amplitudes are usually assumed to 
be negligible as they are helicity suppressed or, equivalently, the form factors are suppressed at 
large q 2 = m 2 D . At first glance, it appears that the factorizable weak annihilation amplitudes are 
too small to be consistent with experiments at all. However, in the diagrammatic approach here, 
the topological amplitudes E and A do receive contributions from the tree and color-suppressed 
amplitudes T and C, respectively, via final-state rescattering (see, e.g., Fig. 1 of Ref. [23|). There- 
fore, even if the short-distance weak annihilation vanishes, long-distance weak annihilation can be 
induced via inelastic FSI's. 

For D — > PP decays, it is expected that the long-distance VF-exchange is dominated by resonant 
FSFs as shown in Fig. 1(a) of Ref. 29J owing to the fact that an abundant spectrum of resonances 
is known to exist at energies close to the mass of the charmed meson. It is easy to draw the 
relevant diagrams, but the difficulty is how to calculate them. In principle, one can evaluate 
the final-state rescattering contribution at the hadron level (see, e.g., Fig. 5 of Ref. |29j), but it 
involves many theoretical uncertainties. As stressed in the Introduction section, a great merit 
of the diagrammatic approach is that the magnitude and phase of the topological TV-exchange 
and VF-annihilation amplitudes can be extracted from the data [see Eq. ([3])]. Moreover, the large 
magnitude and phase of weak annihilation can be quantitatively and qualitatively understood in 



the following manner. As emphasized in Refs. |33l. 1341] . most of the properties of resonances follow 



from unitarity alone, without re gar d to the dynamical mechanism that produces the resonances. 
Consequently, as shown in Refs. [33j, [35(, the effect of resonance-induced FSI's can be described in 
a model-independent manner in terms of the mass and width of the nearby resonances. It is found 
that the E and A amplitudes are modified by resonant FSI's as (see, e.g., Ref. [351 ]) 

E = e + (e 2 * - 1) (e + |) , A = a + (a 2 * - 1) (a + |) , (9) 

with 

e 2i5r = l-i — , (10) 

rriD — fRR + u/2 

where the VF-exchange amplitude E and IF-annihilation amplitude A before the resonant FSI's 
are denoted by e and a, respectively. Therefore, even if the short-distance weak annihilation is 
turned off, a long-distance VF-exchange (VF-annihilation) contribution can still be induced from 
the tree amplitude T (C) via FSI rescattering in resonance formation. To see the importance of 
resonant FSI's, consider the scalar resonance -fTo(1950) with a mass 1945 ± 10 ± 20 MeV and a 
width 201 ± 34 ± 79 MeV, contributing to the IF-exchange in D — > Kir, Krj. Assuming e = in 
Eq. ([9]), we obtain E = 1.68 x 10 -6 exp(zl43°) GeV, which is close to the "experimental" value 
of E given in Eq. (j3|). This suggests that weak annihilation topologies in D — > PP decays are 
dominated by nearby resonances via final-state rescattering. Contrary to the PP sector, we have 
shown in Ref. [29J that weak annihilation in VP systems is dominated by final-state rescattering 
via quark exchange. 
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III. TOPOLOGICAL AMPLITUDES AND QCD FACTORIZATION 



Although the topological tree amplitudes T, C, E and A for hadronic D decays can be extracted 
from the data, we still need information on penguin amplitudes (QCD penguin, penguin annihila- 
tion, etc.) in order to estimate CP violation in the SCS decays. To calculate the penguin effect, 
we start from the short-distance effective Hamiltonian 



H 



eff 



Gp 
V2 



A p (ciOf + c 2 Of + c 8g 8g ) - X b Y CiOi 

p=d,s 



i=3 



where A r 



V* p V up for p 



d,s,b, and 



0\ = (pc) v _ A (up) v _ A , 
03(5) = {uc) V -A^2(qq) 



CP 



a c 0J V -A 



a) V -Ai 



V^Al 



04(6) = {u a C(3) v _ A ^2{qpq a ) 



aJv^A> 



(11) 



°8g = --^mcucr^il + l^G^c , 
with O^-Oq being the QCD penguin operators and {qiq2) v±A 



(12) 



917^(1 ± 75)92- The electroweak 
penguin operators are not included in the Hamiltonian as they can be neglected in practice. For 
the Wilson coefficients, we take c\ = 1.21, c 2 = —0.41, C3 = 0.02, C4 = —0.04, C5 = 0.01, cq = —0.05 
and cg g = —0.06 from Ref. 0| evaluated at the scale [i = m c . 

To evaluate the hadronic matrix elements of the 4-quark operators, it is necessary to specify 
a suitable framework for the task. We will work with the QCDF approach [la . |36j |. As shown in 



26J, various topological amplitudes of D — > M\M 2 decays can be expressed in terms 



details in Ref. 

of the quantities calculated in the framework of QCDF as follows: 

m Gf 



C 



E 



A 



PE 



PA 



—j= a\ (MiM 2 )X^ DMl ,M2 " 1 , 
v2 

^a 2 (MtM 2 ) X (° M ^ , 
v2 

G F 

{ifDfNhfhh) [bl] Ml M 2 ' 

G F 

—^{ifDfMjM 2 )[fo]MiM2 ' 

% [al(M 1 M 2 ) + V r^a p 6 (M 1 M 2 )}X^ DM ^ , 
a/2 

G F 

—^{ifDfM 1 fM 2 )[h] Ml M 2 > 

G F 

—^{ifDfM 1 fM 2 )[h] MlM2 , 



(13) 



where 77 = 1 for M\M 2 
the quantities 6j with i 
Ref. 



= PP,PV and 77 = — 1 for M\M 2 = VP. In the above equations, 
1, • • • ,4 expressed in terms of annihilation amplitudes are defined in 



361 ] . and A is a factorizable matrix element given by 

x (dp u p 2) _ ( p 2 |jM| )(P 1 |j;|L>)=^ P2 ( m 



X (DP,V) = (F|JM| )(P|J'|D) 

x (dv,p) _ (p|jM| )(v|^|D) 



2 1 
m Pt, 

DP 1 



F DP ^mp 2 ] 



2f V m DPc F^(m z v ) , 
2f P m DPc A° v (mp) , 



(14) 
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with p c being the center-of-mass momentum of either final state particle. Here we have followed the 
conventional Bauer-Stech-Wirbel definition for form factors i^oi an d Aq V 37j. For annihilation 
amplitudes, we choose the convention that M\ (M 2 ) contains an antiquark (a quark) from the weak 
vertex. The chiral factors r^ 2 in Eq. (|13j) are given by 



'(/*) 



2mi 



2rn v /y(/x) 



(15) 



m c (/i)(m 2 + mi)(/i)' ' x Vf ~ y m c (n) f v 

with fy(n) being the scale-dependent transverse decay constant of the vector meson V. The flavor 
operators aP i in Eq. ()13p are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions |36l | 



of(MiM 2 ) 



Cj±l 



Ni{M 2 ) + 



<H±i Cfol s 



N c Air 



An 



V i (M 2 ) + —H i (M 1 M 2 ) + Vf(M 2 ), (16) 



N, 



where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), c, are the Wilson coef- 
ficients, C F = (iV c 2 - l)/(2iV c ) with N c = 3, M 2 is the emitted meson, and Mi shares the same 
spectator quark as the D meson. The quantities Vi(M 2 ) account for vertex corrections, Hi{M\M 2 ) 
for hard spectator interactions with a hard gluon exchange between the emitted meson and the 
spectator quark of the D meson, and Vi(M 2 ) for penguin contractions. The explicit expressions of 
Vi, Hi and V% can be found in Ref. [36|]. The quantities Ni{M 2 ) vanish when % = 6, 8 and M 2 = V, 
and are equal to unity otherwise. 

In general, the decay amplitude is evaluated at the scale \x = m c (m c ). However, as stressed 
in Ref. [ijj, the hard spectator and annihilation contributions should be evaluated at the hard- 
collinear scale fj,h = VfJ^h with 500 MeV. This means Hh ~ 0.8 GeV for D — > M\M 2 decays, 
which is beyond the regime where perturbative QCD is applicable. Therefore, we shall not consider 
the spectator contributions to aj. 

Let us first consider the penguin amplitude 



(17) 



with 



al(7rir) 
a^vrvr) 



C4 + 



C6 + 



C3^ 






-K + 


4vr 2 


N C J 


47T 


N c 


C5^ 






-K + 


4vr 2 


N C J 


Ait 


~N~c 



(18) 



The strong phase of the QCD penguin amplitude arises from vertex corrections and penguin con- 
tractions. The order a s corrections from penguin contraction read [36 1 

CFOi s 




4 , m c 2 _ . . 

3 h 7 + r^ w 



+ C 3 



8 m c 4 
-In— + - 
3 [i 3 



Gm 2 i s i 



Gm 2 (1) 



16 m c 
— in — 
3 /i 



Gm 2 ( s u) — GM 2 ( s d) 



Gm 2 (s s 



4 , m c 2 A . , 
-ln- + --G M2 ( Sp ) 



+ c 3 



8 in ■ 4 

3 ln— + 3 ~ Gm 2 {s u ) - Gm 2 {1) 



+(c 4 + c 6 ) 



16, m c 
— In— 
3 /i 



Gm 2 {s u ) — Gm 2 (scl) — Gm 2 (s s ) — GW 2 (1) 



(19) 
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where Sj = mf/m^, 

Gm 2 (s)= [ dxG(s, 1 -x)$m 2 (x), G M2 (s)= [ dxG(s, 1 - x)$ m2 (x), (20) 
Jo Jo 

and G(s,x) = —4 Jq 1 duu(l — u)ln[s — u{l — u)x\. Here <3>m 2 (^m 2 ) is the twist-2 (-3) light-cone 
distribution amplitude for the meson M 2 . 



A. Weak points of QCD-inspired approaches 



Although QCD-inspired approaches such as pQCD [38| and QCDF |39|-|43j have been applied 
to charm decays, their results cannot be taken seriously. Since the charm quark mass is not heavy 
enough, the l/m c power corrections are so large that a sensible heavy quark expansion is not 
allowed. The large magnitude and phase of weak annihilation (E or ^4) topological amplitude, 
which is of order l/m c in heavy quark limit, is an indicative of the importance of power corrections 
in charm decays. For the parameters a\ and a 2 in D° — > K~ir + decays, we find 

ai (Kir) = lA9e i26 ° , a 2 (Kvr) = 0.37 e^ 155 " , (21) 



where the strong phases arise from the vertex correction. In Ref. [29(, we have shown that for 
D — > Kn decays 

|ai| = 1.22 ±0.02 , \a 2 \ = 0.82 ±0.02 , a 2 /a 1 = (0.67 ± 0.02)e^ (152±1) ° , (22) 

as obtained from the topological amplitudes T and C given in Eq. @. It is clear that, while the 
calculated a\ is consistent with "experiment", the predicted a 2 in QCDF is too small in magnitude 
and this again calls for large l/m c corrections. 

While there is no trustworthy theoretical framework for describing the hadronic D decays, in 
this work we shall use the NLO predictions of QCDF to give a crude estimate of the short-distance 
penguin contributions. 

IV. DIRECT CP VIOLATION 
A. Tree-level CP violation 

Direct CP asymmetry in hadronic charm decays defined by 

T(D -> f) - Y(D -> j 
T(D^f)+T(D^f) 



can occur even at the tree level 44|. Take the — > K°ir + mode as an example. Its decay 
amplitude reads X d (T + P d + PE) + \ S (A + P s + PE) (see Table [T]). The interference between the 
color-allowed tree and VF-annihilation amplitudes leads to the tree-level CP asymmetry 



sin 5 at , (24) 



a ^( D + _> K \+) - 2Im ^ A *) lm{T * A) ~ 1 2 x 10- 3 - 
a diri D s^Kn)- |Ad|2 |r _ A|2 -1.2x10 T 

where 5 at is the strong phase of A relative to T and we have taken into account the fact that the 
magnitude of A is much smaller than T. It is obvious that direct CP violation in charm decays is 
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CKM suppressed by a factor of 10 -3 . Assuming the same topological amplitudes in SCS and CF 
decays, we then find from Eqs. ^ and (|24p that a^ ee ^ ~ 10 -4 for P>+ — > K°tt + . Larger direct 
CP asymmetry can be achieved in those decay modes with interference between T and C or C and 
E. For example, is of order —1.1 x 10~ 3 for D° — > urj' and of order (0.6 ~ 0.7) x 10~ 3 for 

D° -> tt 7] and D° -> K°K*°. 



B. Penguin-induced CP violation 



Direct CP violation does not occur at the tree level in some of the SCS decays, such as 
D — > K + K~ and D — > 7r + 7r _ . In these two decays, the CP asymmetry can only arise from 
the interference between tree and penguin amplitudes 

2Im(A d A*) Im[(T* + E* + P d * + PE + PA)(P S + PE + PA)]^ 



-,( lo °p). 

l dir 1 



7T 1 7T 



1.2 x 10" 



IT -i- E I 2 

I J- 7T7T ' ^WIT | 



sin (5„ 



where S n7T is the strong phase of P* 



(loop) 
l dir 



(K + K~ 



T + E 

+ PE n7T + PA wn relative to T nn + E nn . Likewise, 
P d + PE + PA 



-1.2 x 10 



-3 



T + P 



sin^ft: . 



(25) 



(26) 



KK 




To estimate their CP asymmetries, we first discuss the rates of D [> — > it + it and D° — > K + K . 
Experimentally, the ratio R = T(D° -)• K + K^)/T(D° -)• vr+vr") is about 2.8 0, while it should 
be unity in the SU(3) limit. This is a long-standing puzzle since SU(3) symmetry is expected 
to be broken at the level of 30% only. Without the inclusion of SU(3) breaking effects in the 
topological amplitudes, we see from Table [J that the predicted rate of K + K~ is even smaller 
than that of tt + tt~ due to less phase space available to the former. In the factorization approach, 
SU(3)-breaking effects in the tree amplitudes T are given by 



fx F DK (ml 



K) 



U F DK (ml) 



7T7T 



m 2 D 



rat 



\mt 



m 



D 



m\ F DK (ml) ■ 



(27) 



Using the form factor q 2 dependence and input parameters given in Ref. 29(, we obtain 



T KK /T= 1.275 



T n7T /T = 0.96 



This leads to B{D° -»• K + K~) = (3.4 ± 0.1) x 10~ 3 and J3{D° ->■ vr+vr- 



(28) 

(2.1 ± 0.1) x 10" 3 



assuming no SU(3) symmetry breaking in E, i.e., E KK = E nn = E. Therefore, SU(3) breaking in 
spectator amplitudes leads to R = 1.6 and is still not sufficient to explain the observed value of 
R ~ 2.8. This calls for the consideration of flavor symmetry SU(3) violation in the H^-exchange 
amplitudes. We have argued in Ref. [29j] that the long-distance resonant contribution through the 
nearby state /o(1710) could account for SU(3)-breaking effects in the IF-exchange topology. This 
has to do with the dominance of the scalar glueball content of /o(1710) and the chiral suppression 
effect in the ratio r(/o(1710) — > 7T7f)/r(/o(1710) —> KK). To fit the measured branching fractions 
we find 



E, 



1.6 x 10~ 6 e il08 °GeV 



1.05 e 



-«14 c 



E. 



Ptt 



1.3 x 10" 6 e il45 °GeV = 0.85e i23 °P. (29) 
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TABLE III: Predictions of direct CP asymmetries in units of 10 -3 , where o!^^ denotes CP asym- 
metry arising from tree amplitudes. The superscript (t + p) denotes tree plus QCD penguin am- 
plitudes, (t + pa) for tree plus weak penguin annihilation (PE and PA) amplitudes and "tot" for 
the total amplitude. For the VP modes, we use the solutions (A) and (Al) given in Eq. (jHJ) for tree 
and weak annihilation amplitudes. For weak penguin annihilation, we assume that it is similar 
to the topological E amplitude [see Eq. pip]. World averages of experimental measurements are 
taken from Ref. 0|. Due to the lack of information on the topological amplitudes Ap and Ay, no 
prediction is attempted for D + — > VP and Df — > VP decays. 



Decay Mode 


(tree) 


(t+ P ) 

u dir 


(t+pa) 
a dir 


_(tot) 
u dir 


Expt(10- 3 ) 


Decay Mode 


(tree) 


(t+P) 

a dir 


(t+pa) 
a dir 


(tot) 
a dir 


D a - 


■> 7T+7T 





-0.04 


0.90 


0.86 


2.0 ±2.2 


D°- 


■> 7T + p 





0.10 


-0.60 


-0.51 


D°- 


> 7r°7r° 





0.25 


0.59 


0.85 


1±48 


D°- 


-> ir~ p + 





-0.05 


-0.22 


-0.27 


D°- 


> 7T°?7 


0.63 


0.43 


0.03 


-0.16 




D°- 










-0.74 


-0.74 


D° - 


>7rV 


-0.51 


-0.60 


0.09 


-0.01 




D°- 


> K+K*- 





-0.10 


0.60 


0.50 


D a - 


* w 


-0.37 


-0.38 


-0.70 


-0.71 




D°- 


* K-K*+ 





0.06 


0.22 


0.29 


D° - 




0.39 


0.44 


0.21 


0.25 




D°- 


■> K°K*° 


0.73 


0.73 


0.73 


0.73 


D°- 


> K+K- 





0.02 


-0.50 


-0.48 


-2.3 ± 1.7 


D°- 


■* K°K*° 


-0.73 


-0.73 


-0.73 


-0.73 


D+ - 


-> 7T+7T 
















D° - 


+ n°u 





-0.15 


0.53 


0.37 


D+ - 


-> TT + r/ 


0.37 


0.28 


-0.56 


-0.65 


17.4 ± 11.5 a 


D° - 


■> 7T> 














D+ - 


+ 7r + r/ 


-0.21 


-0.26 


0.42 


0.41 


-1.2 ±11.3 a 


D° - 


■> 


0.19 


0.19 


0.50 


0.50 


D+ - 


-> K+K° 


-0.07 


0.08 


-0.53 


-0.38 


-0.9 ±6.3 


D°- 


■> Tj'ui 


-1.07 


-1.05 


-0.91 


-0.89 


Df- 




0.09 


-0.07 


0.69 


0.52 


65.3 ±24.6 


D° - 


-> r](j) 














Df- 




0.01 


0.02 


0.87 


0.88 


20 ± 290 


D°- 




-0.53 


-0.54 


-0.22 


-0.23 


Df- 




-0.61 


-0.27 


-0.53 


-0.19 


-200 ± 180 


D°- 


+ 77V 


0.59 


0.58 


0.21 


0.20 


Df~ 




0.35 


0.58 


-0.63 


-0.41 


-170 ±370 















a Data from 48 



Using the input parameters for the light-cone distribution amplitudes of li ght mesons, quark 



masses and decay constants from Refs. 45j, |46|] and form factors from Refs. (29|, |47j], 2 we find 
ps \ / pd \ 



T + E 



0.35e- il75 °, 



T + E 



0.24 e 



-il76° 



(30) 



KK 



Hence, 5-n 



t+p 

l dir 



(K + K~ 



S KK = -176°. From Eqs. QEft and ((26]), we derive ^^(vr+vr") = -3.9 x 10~ 5 and 
= 2.0 x 10 -5 , where denotes the CP asymmetry arising from the interference 
between tree and QCD-penguin amplitudes. QCD-penguin induced CP asymmetries in D° — > 
7r + 7r _ , K + K~ are small partially due to the almost trivial strong phases of 8^ and 5kk- 

Short-distance weak penguin annihilation contributions (PE and PA) to charm decays are very 
small; typically, PE/T and PA/T are of order 10 -5 and hence are negligible. Nevertheless, long- 
distance contributions to SCS decays, for example, D° — > it + tt~, can proceed through the weak 
decay D° — > K + K~ followed by a resonant-like final-state rescattering as depicted in Fig. [2 It 



2 More specifically, we use fjj = 219 MeV and fxj s = 260 MeV [46j for the charmed meson decay constants 



and those in 



47] for D — > V transition form factors. 
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A 



D< 




FIG. 2: Long-distance contribution to D° — > ir + ir through a resonant-like final-state rescattering. 
It has the same topology as the QCD-penguin exchange topological diagram PE. 



has the same topology as the QCD-penguin exchange topological graph PE. Just as the weak 
annihilation topologies E and A, it is expected that weak penguin annihilation will receive sizable 
long-distance contributions from final-state interactions as well. Recall that soft corrections due 
to penguin annihilation have been proposed to resolve some problems in hadronic B decays, for 
example, the rate deficit problem for penguin-dominated decays and the CP puzzle for B° — > K~ir + 



361 ] . Hence, we shall assume that PE, PEp and PEy are of the same order of magnitude as E, Ep 



and Ey, respectively. For concreteness, we take (in units of 10 6 ) 

PE = 1.6e ill5 °GeV, PE P = 3.1 e~ m ° , PEy = 1.5e" il20 °. (31) 



C. Numerical results and discussion 

The calculated direct CP asymmetries for various SCS D — > PP and D — > VP decays are 
summarized in Table IIIII 3 Before embarking on CP violation, we have taken into account the 
SU(3) breaking effects in tree amplitudes T and C using the factorization approach so that the 
predicted branching fractions are consistent with the data shown in Tables U and HH Besides the 
aforementioned example for D° — > tt + tt~ and D° — > K + K~, we found, for example, T n+ „o/T = 
0.85, „o/C = 0.93 for the decays D + — > 7r + ?/') 29j] . The predicted branching fractions 



B(D+ tt+7]) = (3.02±0.19)xl0~ 3 and B(D+ tt+t]') = (4.69±0.21) x 10~ 3 with A n+ (>) /A = 0.7 
are in agreement with the measured ones B{D + — > tt + t]) = (3.53 ± 0.21) x 10 -3 and B(D + — > 
vr+r/) = (4.67 ±0.29) x 10" 3 (see Table U). 

In general, tree-level CP violation lies in the range 10" 4 < o£ ee) < 10~ 3 . The largest one 
occurs in the decay D° — > ujt( where a^j, = —1.1 x 10~ 3 . We would like to accentuate once 



3 



CP -violating asymmetries in SCS decays of charmed mesons have also been investigated in Ref. [49 1. 



In this work, final-state interaction effects are studied by assuming that FSFs are dominated by nearby 



resonances, similar to the work of Refs. |33l. |35| . The major uncertainties there arise from the unknown 
masses and widths of the resonances at energies near the charmed meson mass. Noticeable differences 
between Ref. 49( and our results for direct CP asymmetries are that Aa^'p of D° — > K* + K~ and 
D° — > K*~K + arc predicted to have the same sign in Ref. 49( and likewise for D — s- K + K~ and 
D° — > 7r + 7r~, whereas they are of opposite signs in our work. 
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again that the estimation of tree-level CP asymmetries in the diagrammatic approach is reliable 
and trustworthy since the magnitude and phase of the tree topological amplitudes are extracted 
from the measured data. By an inspection of a£ e , one can get an idea about the size of direct 
CP violation within the SM. 

The QCD penguin effects on CP asymmetry can be inferred from a comparison between a^ p ^ 
and a^ e \ We see from Table Hill that the short-distance penguin effects are small compared to the 
tree-amplitude-induced CP violation. For decay modes with vanishing o£ , the QCD-penguin- 



induced CP asymmetry is of order (1 ~ 2) x 10~ 4 or even smaller (see Table ITTTj) . 

When the short-distance weak penguin annihilation contribution is estimated in QCDF at the 
hard-collinear scale Hh ~ 800 MeV, it is found to be very small as mentioned before. But weak 
penguin annihilation has the chance to be greatly enhanced. In QCDF, it involves endpoint diver- 
gences related to the soft gluon effects. In the diagrammatic approach, it receives long-distance 
contributions from the color-allowed tree amplitude followed by final-state rescattering as shown 
in Fig. [2j Assuming that penguin annihilation is similar to the topological VF-exchange amplitude 
[see Eq. (|3"T]) ]. its effect can be studied by comparing a^ pa ^ with a^ r , where the superscript 
"pa" denotes weak penguin annihilation. It is evident from Table IIIII that penguin annihilation 
could play an essential role in the study of CP violation. 

After summing over all possible contributions, we see that the predicted CP violation denoted 
by (or a^p) is at most of order 10 -3 in the SM. For Aa^p, the CP asymmetry difference 
in D° — > K + K~ and D° — > 7r + 7r _ , we obtain a value of —0.13% which is too small compared 
to the world average of —(0.645 ± 0.180)%. Since in the SM, Aa^Jp stems mainly from weak 
penguin annihilation, we can vary the amplitude PE to see how much enhancement we can gain. 
Even with the maximal magnitude \PE\ ~ T and a maximal strong phase relative to T, we get 
Aa$? P = -0.25%. This is more than 2a away from the current world average. Hence, if the LHCb 
result for Aa^Jp is confirmed by further data analysis, it will imply new physics in the charm sector. 



V. CONCLUSIONS 



In view of the recently measured CP asymmetry difference Aa^Jp between the D° — > K + K~ and 
D° — > 7r + 7r~ modes, we scrutinize the direct CP violation in the singly Cabibbo-suppressed (SCS) 
D — > PP and VP decays within the standard model (SM), where P and V refer to pseudoscalar 
and vector mesons, respectively. Such an analysis is helpful in diagnosing possible evidence of new 
physics in the charm sector. 

Direct CP violation in such D decays may arise from the interference between tree- 
level amplitudes of different topologies. However, the magnitude is proportional to 
lm{V* d V ud V* s V us )/\V* d V ud \ 2 , which is O(10~ 3 ). Modulated by the size ratio of different flavor 
amplitudes and the sine of the relative strong phase, we expect that such tree-level CP asymmetry 
should be at the order of 10 -3 or even smaller. The great merit of the topological approach is 
that the magnitude and the relative strong phase of each individual topological tree amplitude in 
charm decays can be extracted from the data. Hence, the estimate of a^ p should be trustworthy 
and reliable. We find that a^p(D° — > wq') ~ —1.1 x 10~ 3 is the largest one among all the SCS 



1(3 



modes. 

In the D° — > K + K~ and D° — > tt + tt~ decays, the tree-level amplitudes have exactly the same 
CKM factors and thus do not induce direct CP asymmetry. In such cases, one must invoke the 
penguin amplitudes, including the QCD penguin and weak penguin annihilation amplitudes. In 
addition, we have taken into account the SU(3) breaking effects in the color-allowed tree T diagram 
and the resonance effects in W-exchange E diagram, as done in Ref. [29[], in order to explain the 
observed branching fractions of the two modes. As given in the a^ p ^ column of Table 1111} the 
direct CP asymmetries of both modes are at a few x 10~ 5 level. This is seen to be largely due to the 
trivial relative strong phase between the QCD penguin amplitude and the tree-level amplitudes. 

We next include the contributions of weak penguin annihilation diagrams, PE, PEp and PEy . 
The short-distance contributions of such diagrams are typically five orders of magnitude smaller 
than the color-allowed tree diagram and thus negligible. However, they may be enhanced by 
long-distance final-state rescattering through resonances, as in the case of the weak annihilation 
diagrams E and A. We manage to maximize the CP -violating effects from the interference with 
these weak penguin annihilation diagrams by assuming these diagrams are as large as the tree-level 
E, Ep, and Ey, respectively, and have sizable strong phases relative to T. In the end, we observe 
that the magnitude of CP asymmetry in the SCS charm decay modes can reach at most 10~ 3 . In 
particular, the CP asymmetry difference Aa^p between D° — > K + K~ and D° — > 7t + tt~ is found 
to be about —0.13%. Even with the parameter choices of \PE\ ~ |T| and a relative strong phase 
of 90°, we find Aa^p = -0.25%, which can be regarded as the upper bound on Aa$ P in the SM, 
still more than 2a away from the world-average experimental result —(0.645 ± 0.180)%. 

In conclusion, if Aa^Jp continues to be large with more statistics in the future or if the direct 
CP asymmetry of any of the discussed modes is significantly larger than 10~ 3 , it will be clear 
evidence of physics beyond the SM in the charm sector. 
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